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We propose an effective model of strongly coupled gauge theory at finite tem- 
perature on R 3 in the presence of an infrared cutoff. It is constructed by 
considering the theory on S 3 with an infrared cutoff and then taking the size 
of the S 3 to infinity while keeping the cutoff fixed. This model reproduces 
various qualitative features expected from its gravity dual. 
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1 Introduction 



Considerable progress has been made recently in understanding various aspects of gauge 
theories with confinement [1,2,5-9]. It turned out that one way of fabricating gravity dual 
of a gauge theory with confinement is to introduce an ultraviolet cut-off. Indeed, in [10], 
the gravity dual of the confining gauge theory on flat space was constructed by smoothly 
capping off the geometry at small radius. A simpler model appeared in [12,13] where the 
gravitational description consists of AdS 5 with small radius region removed and this was 
termed as a hard wall approximation. It is shown [13] that for small enough cutoff (infact, 
smaller than the size of a black hole horizon), one obtains Hawking- Page transition at a 
finite nonzero temperature. Introduction the cut-off in a milder way leads to a soft wall 
approximation where one uses an appropriate dilaton field in the background to truncate 
small radius region. However, in either cases, the short distance cutoff on the gravity side 
translates into an IR cutoff in the corresponding gauge theory and so the gauge theory 
dual to the gravity theory under consideration is a theory on R 3 in the presence of an IR 
cutoff. Using the gravity dual one can learn various aspects of this gauge theory in the 
strongly coupled regime. At this point one may ask: Can these lessons be used to design 
an effective theory for the strongly coupled regime of this gauge theory with IR cut-off? 
This is essentially the point that we would like to address in the present note. 

Necessary hint comes from the consideration of the confining gauge theory on a 
compact space, such as, S* 3 . There one does not need a cut-off and radius of the three- 
sphere provides the scale essential for confinement. Indeed an effective model for strongly 
coupled regime of this gauge theory on S 3 was constructed in [2,4]. This model, which 
we call (a, b) model in the sequel, is characterized by two parameters a and b. These 
parameters, in general, depend on temperature and the 't Hooft coupling. If a, b satisfy 
certain properties, this model reproduces the qualitative features of the phase structures 
expected from the dual theory in the bulk. To construct an analogus model on R 3 with 
an IR cutoff, we use the following route: 

• Incorporate the IR cut-off in the (a, b) model so as to describe a gauge theory 
on S* 3 with an infrared cutoff. 

• Keep the cutoff fixed and let the size of S 3 go to infinity. 

The first step of constructing an intermediate model is easier. On generic ground, we 
expect that, in the presence of cutoff, the parameters a and b would depend on the cutoff. 
Indeed, we find that the (a, b) model reproduces the qualitative features of gauge theory 
on S 3 with cutoff if a, b depend on the cutoff in a specific way. This can be computed 
by comparing the effective potential in the gravity side with the free energy of the gauge 
theory arising from the (a, b) model. In the absence of cutoff, a and b are monotonically 
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increasing function of the temperature [2] . However, we find, while b remains an increasing 
fuction even in the presence of cutoff, the nature of the variation of a is sensitive to the 
magnitude of the cutoff. The second step is more difficult. Since the radius of the three- 
sphere does not appear as a tunable parameter, there is no simple procedure to carry out 
this step. Comparing gravity duals of the gauge theories on S* 3 and R 3 respecvtively we 
propose that the effective model remains qualitatively similar. Based on this assumption 
we have constructed an effective model. Details of the computation is presented in latter 
sections of the note. Indeed our model posseses all the qualitative features expected 
from the gravity side. However, it also leads to a puzzle which we address in the next 
paragraph. 

One of the features of our effective model is that there always exists an unstable 
maximum which separates the two minima. These two minima, in the gravity side, 
represent thermal AdS and black hole in presence of UV cutoff. On the other hand, one 
expects the maximum to represent a small unstable black hole. However, for the cases 
analysed in [13], such a black hole is missing. We suspect that this is due to the fact that 
the cutoff was introduced by hand. Once this cutoff is intruduced, both in soft and hard 
wall models dynamically, such black hole would prehaps appear. A preliminary support 
of our proposal comes from a calculation carried out in the apendix of this paper. Here we 
analyse a gravity action which has a fine tuned dilaton potential. If we set the dilaton to 
zero, we have a flat AdS black hole. The temperature of this black hole increases linearly 
with the horizon size. However, as we turn on a dilaton profile and the dilaton potential ( 
satisfying the BF bound), we see that for a given temperature, there can have two black 
holes - one of which is unstable. 

The plan of the paper is as follows. In section (2), we analyze the gravity theory 
with a radial cutoff. In section (3), we consider the zero coupling limit of gauge theory 
with an IR cutoff. In section (4), we study the behaviour of the phenomenological matrix 
model on S 3 with an IR cutoff. Next section, (5), consists of analysis of the model as we 
take the size of S* 3 to infinity, keeping the cutoff fixed. This paper ends with a discussion 
of our results. Finally, in the appendix, we show, with in the soft wall model, how an 
unstable black hole may appear if we introduce cutoff dynamically. This is a prediction 
that arise from an analysis of our martix model of R 3 . 

2 Bulk with UV cutoff 

We begin this section, with gravitational configuration having boundary S 3 x S* 1 and its 
thermodynamic properties the presence of an UV cutoff. In the latter part of this section, 
we take a suitable limit to obtain gravitational configuration having boundary R 3 x S* 1 
and study its analogous properties. 
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• Gravitational configuration with boundary S 3 x S* 1 



Five dimensional gravity action with a negative cosmological constant is given by 

where k 2 is the five dimensional gravitational constant and I 2 parametrizes the cosmolog- 
ical constant. This action admits two solutions with same asymptotic behaviour: AdS 
and AdS black hole. For both, the metric can be written as 

dS 2 = -V(r)dt 2 + V(r-y l dr 2 + r 2 dQ 2 3 . (2.2) 

While for AdS, V(r) is 



V(r) = 1 + (2.3) 



for the AdS black hole it has the form 



r 2 m 



V(r) = 1 + ^ - - 2 . (2.4) 

In the last equation m is associated to the ADM mass of the black hole and is given by 

M = (2.5) 

where u>3 is the volume of the unit three-sphere. In the asymptotic limit, r — > oo, the 
metric becomes the AdS metric. The black hole has a curvature singularity at r = and 
has a single horizon. The horizon radius r + can be obtained by setting V(r) = and is 
given by the real positive root of 

l-£ + £ = 0. (2.6) 

Various thermodynamic quantities for the black hole with an UV cutoff is given 
in [14]. Continuing to Euclidean signature, with time direction as a circle, one finds that 
smoothness of the geometry at non-zero r requires the radius of the time circle to vanish 
at horizon. This leads to a temperature, 

2r 2 4- 1 2 

which is independent of the cutoff. There is a critical temperature To, termed as nucleation 
temperature, below which there is no black hole. For T > Tq there are two black hole 
solutions - an unstable small black hole and a stable big black hole with horizon radius 
given by the smaller and larger root of r + respectively. They merge at T = T . 



3 



The free energy of the system is obtained by integrating the regularized action for 
AdS space and the black hole over the amputated space where a small three-sphere of 
radius ro is removed from the origin and so the lower limit of the integration on the radial 
coordinate is r = r instead of the origin r = 0. Imposing same asymptotic boundary 
condition for both the geometries and choosing the free energy of the AdS to be zero, one 
can obtain the free energy associated with the black hole as [14] , 

F =-2^( 2r -- r *- 2r »- r +' 2 )' (2 ' 8) 

Here we choose 

r max = r+ for r < r+ (2.9) 

r m ax = r for r > r+. (2.10) 

Hawking-Page (HP) transition occurs only if r < r + as only in this case the free energy 
(2.8) changes sign. In terms of Wilson loop operator, one can see that unless r < r+, 
there is no world-sheet belonging to zero homotopy group which ends on the temperature 
circle at the asymptotic infinity. Consequently, Wilson loop operator does not show a 
phase transition. The same conclusion is reported in [13]. At this point it may appear in 
mind that since cut-off is within the horizon thermodynamic quantities are not affected 
by it. Indeed we will see the critical exponents do not depend on the cut-off. 

Since we are interested in comparing the thermodynamic behaviour with that of the 
matrix model counterpart coming from gauge theory, where the order parameter is the 
moment of the vev of the Wilson loop operator, we restrict ourselves to r < r + . The 
situation for r > r + is considered in [14]. The entropy of this sytem turns out to be 

S = ^± (2.11) 

and does not depend on the cutoff. 

In order to summarize different phases in a convenient way ( which will help us 
in later sections to make a comparison with the matrix model ), we construct a Landau 
function. This reproduces equations (2.7, 2.8, 2.11) on shell. First, we define the following 
dimensionless quantitites 

r -f = f ) r -f=f^Tl = T^ = R\ (2.12) 

Then, by identifying the horizon radius f as the order parameter, it is easy to write down 
the Landau function J-{f, fo,T) as, 

F{r, f , f) = H (V - Anf f 3 + 3f 2 + 2fg) , (2.13) 
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where jF(f, f ,T) = lJ-(f, fo,T). On shell, this function reduces to the expression of free 
energy given in (2.8). In order to see this, we need to compute the on shell temperature, 
which can be obtained by setting first variation of Landau function with respect to the 
order parameter f equal to zero. The result matches with (2.7). Then substituting it in 
(2.13) we get (2.8). We plot the Landau function for various temperatures with three 
different values of cutoff. The critical temperature corresponds to having two degenerate 
minima of T . This can be determined as follows. The function T vanishes when the 
right hand side in equation (2.13) is identically equal to zero. This, in turn, gives an 
equation for T involving f and fo. Equating this with (2.7), we get f in terms of fo. Upon 
substituting this f in (2.7), we finally get the critical temperature as 




Figure 1: The plot of the Landau function T with r for various temperatures with ro = 0.2, 0.3 
and 0.4 from left to right. 



From the plots we observe that the transition temperature increases with the cutoff. 
We will see that a similar behaviour occurs on the gauge theory side as well. Neverthe- 
less, the qualitative features remain independent of the cutoff. In addition, we find that 
the critical behaviour around the nucleation temperature (To), is universal and does not 
depend on details of the model. A computation of the specific heat near the nucleation 
temperature, T shows that it diverges as (T — Tq)^ 1 / 2 . The critical exponent is thus 
unaffected by the presence of the cutoff. 

• Gravitational configuration with boundary R 3 x S 1 

Having discussed the gravitational configurations with boundary S 3 x S 1 in the 
presence of an ultraviolate cutoff, we now review the case where the boundary is R 3 x S 1 . 
This can simply be obtained by appropriate scaling of the coordinates and the parameters 
that are involved in the previous geometries. We scale 

I 2 A 4 r 4 
r -> Ar, t -> -t, r -> Ar , m -> — K (2.15) 
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Figure 2: The plot of the Landau function T with r for various temperatures with r = 0.1. 



Then, taking the limit (A//) — > oo, earlier black hole metric reduces to 

; (r 2 - T l) de + (r 2 - ^) ~V + r 2 £ <**?] • (2-16) 



»=i 



where Xi are the appropriate coordinates on R 3 . The horizon radius in this limit reduces 
to r h . 

Since the temperature is related to the inverse of the Euclidean time, in this scaling 
limit, the temperature of the black hole reduces to 

T=— . (2.17) 

7T 

We note that, unlike the case of black hole with spherical horizon, here for a given tem- 
perature, there is only a single black hole with a non-zero horizon. The free energy of this 
flat black hole follows from the spherical ones in the scaling limit and is given by 

F = -^(2rL x - 4 - 2r 4 ), (2.18) 

where 

W = r h , for r <r h , 
= r , for r > r h , 

and we have scaled the free energy A 4 F — > F. The Hawking-Page transition when r$ < rh 
and at temperature T c = 2^(r /n). 

Finally, the Landau function can be obtained easily as 

T^(3f 4 -47rff 3 + 2f 4 ), (2.19) 



where we have defined r = r^/l. For r = 0.1 the plot of this function for different 
temperature is shown in figure (2). 

We end this section with the computation of the potential for a quark- ant iquark 
pair on R 3 x S 1 . The confined phase of the gauge theory corresponds to a thermal AdS 
geometry in the bulk. We thus compute the Wilson loop by computing the area of the 
dual world sheet extending into the bulk AdS space following [15, 16]. The computation 
will however be done in the pesence of the cutoff 1 . 

The AdS metric with boundary R 3 x S 1 is, 



ds 2 = a ^(-dt 2 + dx 2 )+R 2 ^ + R 2 dniy (2.20) 
R 2 = ^/ArtgN. (2.21) 



The AdS radius is I 2 = a R 2 . 



To compute the potential between a quark and antiquark pair placed at x — ±L/2, 
we consider a rectangular Wilson loop in the x — t plane. The dual area can be computed 
from the Nambu-Goto action, 



S = — ^ / dadr V 'det (G M Nd a X M d T X N ), (2.22) 
2na J 

where, Gmn are the components of the metric (2.20). In the static gauge: a = x and 
r = t, we get 



S= ^j dx ^MW^WIW- (2-23) 



Equation of motion from (2.23) satisfies, 

U 4 



vW + {U/RY 



C, (2.24) 



where C is a constant that can be fixed by saying that d x U = for U — U (U is the 
cutoff radius). So that we get C = R 2 Uq. 



1 Similar computations, with the bulk AdS Schwarzschild black hole geometry, were carried out in 
[17, 18] where phase transition different from the one studied in the present paper, was considered. 
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Figure 3: World-sheet configuration dual to the Wilson loop in the presence of cutoff, U . 
Equation (2.24) can be solved for x using the boundary condition that x(U = U ) = 



X — Xq 



R 2 f u ^ dy 



r 

= u A 



(2.25) 



where y = U/Uq. For U = oo, x = L/2. So given L and xq, one can solve for Uq from 
(2.25). The full energy of the configuration is given by the contribution from the sections 
of the curve "AB" plus "BC", shown in figure (3). 

The contribution from "AB" is given by inserting the solution (2.25) into (2.23), 



Eab 



Uo 
2ir 



-l)dy-l 



(2.26) 



This contribution goes as the inverse power of L, because from (2.25), U ~ R 2 / (L/2—x ). 
"BC" is the straight vertical line at U — Uq, the contribution from which is given by 
evaluating (2.23) along this line, 



1 U 2 
E BC = ^[L-L ].f 2 , 

where, in equation (2.27), L /2 = L/2 — x . The string tension is thus, 

4tt i? 2 ' 



(2.27) 



(2.28) 
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Going to the original radial coordinate by replacing U — r/a , 

1 r 2 

This shows that the tension grows with the cutoff. For Af = 4 SYM on S* 3 we do not 
expect to get such a linear potential and hence a string tension. This potential has been 
recently computed in [19]. Though the nature of phase transition in the two cases - that 
is one on finite S 3 and the other on R 3 with a cutoff are similar, the behavior of the 
potential before and after the transitions seem to be quite different in the two cases. 

Let us note that the string tension vanishes as the cutoff is taken to zero in equation 
(2.29). This is expected because, in this case, the theory is same as that of one on 
the three-sphere with the radius of the sphere taken to infinity where we only have the 
coulombic term in the potential. Finally, if we consider the theory on S 3 with an additional 
cutoff, this R 3 result can also be thought of arising when L/R S 3 << 1, where -R53 is the 
radius of the three-sphere. 



3 Effective model at zero coupling with IR cutoff 



In this section, we discuss an effective model obtained from gauge theory on a three-sphere 
in presence of an IR cut-off. A three-sphere provides a natural scale which translates into 
confinement scale and in view of that it may appear superfluous to consider an additional 
cut-off. However, as emphasized in the introduction this model is an intermediate step 
to construct an effective theory for strongly coupled gauge theory in R 3 in presence of 
an IR cut-off. Before we move to the discussion of the boundary gauge theory at strong 
coupling, in this section we will begin with the partition function of M — 4 U (N) super 
Yang-Mills theory at zero coupling on a three-sphere with an IR cutoff following [4] . The 
field content of the theory consists of vector bosons, scalar bosons and fermions all in 
adjoint representation of U(N). We expand various fields into spherical harmonic modes 
on the three-sphere leading to a quantum mechanical system consisting of infinite number 
of degrees of freedom. The non-zero modes lead to massive degrees of freedom and it turns 
out that the zero modes are also massive in this problem. The fermionic zero modes are 
massive because of anti-periodic boundary condition around the temperature circle and 
the scalar zero modes are massive due to their coupling with the Ricci scalar of the three- 
sphere. Since we set the coupling to zero these degrees of freedom can be thought of as 
infinite number of decoupled harmonic oscillators each in the adjoint representation of 
U(N). They can be either bosonic or fermionic depending on their parent field. Any 
excitation of this system, then, composed of single particle states corresponding to the 
various oscillator modes and a general physical state can be obtained by acting with an 
arbitrary collection of these oscillators on the Fock space vacuum. However, the Gauss's 
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law does not allow any non-singlet state to be physical. So while considering physical 
excitations we need to consider those which are singlets of U(N). 



We begin with the partition function for bosonic modes in adjoint representation 
of U(N). Single particle modes are denoted by % with energy A generic state is 
represented by a set of integers {n,} where n, is the level of i'th mode and so with energy 



niEi. The corresponding partition function can be written as 

i=0 



zw=nE 



X 



riiEi 



S(ni), 



(3.1) 



i n;=0 



where S{rii) represent number of singlet states in a symmetric product of number of 
particles in adjoint representation, which takes care of the fact that we are summing over 
physical states only. Since number of singlet in any representation can be obtained by 
integrating the character of that representation over the group manifold, we write 



S{n)= J [d[/]Xsym"(ad.j)- 

We substitute this integral expression in (3.1) to obtain 

/oo 
-■ n i= o 



(3.2) 



(3.3) 



As shown in the appendix of [4], this expression reduces to 



Z{x) = J [dU] 



exp 



J2 -x mEi tr(U m )tr((tf) m ) 



m=l 



(3.4) 



where we use Xadj = tr(U)tr(W) and the tr implies trace in the fundamental represen- 
tation. We can carry out a partial sum over all single particle states associated with a 
particular field. There could be different single particle partition functions for scalars and 
vectors. For the time being we assume only one kind and write single particle partition 
function as z(x) = ^x Ei . The above partition function in (3.4)reduces to 



r r 00 1 

Z(x)= [dU]exp V — z B {x m )tr{U m )tr((U ] ) m ) 
J i m 

m=l 



(3.5) 



A similar computation can be done for the fermionic fields also (see [4] for details) where 
the symmetric product in (3.1) is replaced by antisymmetric product. Considering both 
the bosons and the fermions we obtain the total partition function as 



. r oo 

Z(x) = / [dU\ exp —\ z B{x m ) + (-l) m+1 z F (x m )]tr(U m )tr((U^) m ) 

J m=l m 



(3.6) 
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where zb and zp are bosonic and fermionic single particle partition function respectively. 
If we have ns number of bosonic fields and % number of fermionic field the expression in 
the square bracket in (3.6) will be replaced by riBZB(x m ) + (— l) m+1 riFZF{x m ). In addition, 
if there are two kinds of bosonic fields, such as, n v number of vector fields and n s number 
of scalar fields the term becomes n v z v (x m ; R) + n s z s (x m ; R) + (— l) m+1 riFZF(x m ) where z v 
and z s are single particle partition function for vector and scalar fields. 

In presence of a cutoff e we will consider the states, with energy satisfying the bound 

oo 

TiiEi > e . The corresponding partition function can be written as 

i=0 

oo 

Z(x; e ) = J] x^OimEi - eo)S(m), (3.7) 



i n;=0 



where we use the ^-function satisfying 9(x) = for x < and 9(x) = 1 for x > 0. Using 
the integral representation 

oo 

0(x) = J y e*°, (3.8) 



we can rewrite the partition function as 

Z(x; co) = / T e^ ]J ^ (xe^Sim), (3.9) 
J ' i n i= o 
rlk 

^e- lkt0 Z{xe lk ). (3.10) 

K 

Then following the similar procedure as above we can write the partition function with 
cutoff from (3.4) by replacing z B and z F by z^ and Zp respectively 



r r 00 1 

Z(x; eo) = / [dU\ exp ^ —\z^{x m ) + {-l) m+1 z^{x m )}tr(U rn )tr((U) ] ) m ) 

J [m=i m 



, (3.11) 



where z e £ and are bosonic and fermionic single particle partition functions respectively 
with cutoff e . Both of them have the general form 



. i f elk 
z eo ( x ) = / — e 



ikeo z(xe ik ). (3.12) 



A formal path-integral derivation of the partition function for Yang-Mills theory at 
zero coupling is done in [4] and turn out to give rise to the same partition function as 
in (3.6). All the modes once again gets massive and can be systematically integrated 
out. The only massless mode that remains in the spectrum is the zero mode associated 
with A around the thermal cycle. Thus the group theoretic quantity j^Tr{U) actually 
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represents Polyakov loop operator around thermal cycle and serves as an order parameter. 
In presence of cutoff the generalization of computation of partition function in a path- 
integral formulation is very similar and so we are not repeating that here. One will get 
the same result as (3.6). 

The only variable in this effective action is an N x N unitary matrix and thus we 
have arrived at an effective matrix model. Following [4] we rewrite it in terms of the 
eigenvalues of U as {e ai } where ctj's are within the range (—it < oli < n) and obtain a 
system of eigenvalues distributed over a circle. The measure can be rewritten in terms of 
the eigenvalues as 

TT 

[^]-n/[^ri sin2 (^i^)- ( 3 - 13 ) 

— TT J 

We make the substitution 

tr(CT) -> ^e ina \ (3.14) 

3 

in the action and write down the partition function as 




(3.15) 



where 

V eo (a) = -log|sin(a/2)| -^2~[z^(x n ) + (-l) n+1 z e ° (x n )\ cos(na). (3.16) 

n=l n 

The temperature independent repulsive potential is independent of cutoff and the attrac- 
tive potential, which monotonically increases with temperature, decreases as the value 
of cutoff increases. That is what we expect because increasing the cutoff amounts to 
discarding lower powers of a; in a Taylor series expansion of zb(x) and zp(x) and so the 
attractive potential dominates over the repulsive potential at higher value of x. 

This system of eigenvalues shows two different phases in terms of the eigenvalue 
distribution as order parameter. At low temperature the repulsive potential dominates 
and the eigenvalues tend to be distributed uniformly over the circle —n < a < it. As 
temperature increases the attractive potential gets stronger and at some critical temper- 
ature the eigenvalues tend to clump together signalling a phase transition as happens in 
absence of cutoff. Since attractive potential decreases with cutoff while repulsive potential 
is independent of it, the phase transition temperature increases. Introducing the eigen- 

7T 

value distribution function, p(a), normalized as f dap(a) = 1 the effective action can be 
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written as 

C r N 2 °° 

S[p(a)] = N 2 / da / dfip{*)p{p)V*{a-P) = — ^|p n |V^, (3.17) 

where p n := J dap(a) cos(na) represents moments of p and the moments of potential 
V n = f daV(a) cos(na) are 

V*° = — (1 - 4°(x n ) - (-l) n+1 4°(x n )). (3.18) 

From (3.17) we see the uniform distribution is minimum when V*° > 0. As we increase 
the temperature it is V\ which reaches the zero first and so the uniform distribution 
will be destabilized. Setting V^°{xh) = we get at x = xh z b( x h) + z f( x h) = 1- 
This temperature is identified with the Hagedorn temperature of the corresponding gauge 
theory. The coefficients zb and zp depends on the bosonic and fermionic matter contents 
respectively. In particular, for M = 4 SYM, in which we are interested, the single particle 
partition functions are given by 

. , 6x + 12x 2 -2x 3 . . 16a; 3 / 2 

MX) = (1 _ x)3 > *F(X) = ^—y 3 . (3.19) 

From this expression, using (3.1) we can calculate and Zp respectively. We have 
plotted the Hagedorn temperature for free M = 4 SYM against the cutoff and observe 
that the temperature decreases with increase in cutoff. 

X 

1 - 

0.7 - 
0.5 - 

0.3 - 
0.2 - 

1 5 10 50 100 50QL000 n 

Figure 4: The plot of xh = e~^l TH ^ with cutoff eo = n. For convenience we have used a LogLog 
Plot. 



4 Effective model at strong coupling on S* 3 with IR 

In the last section we briefly reviewed the derivation of the matrix model associated 
with an M = 4 SYM theory with an IR cutoff on a three-sphere at zero coupling. In 
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order to compare with gravity we need to go to the strongly coupled regime. Let us 
begin considering the analogous analysis [2,4] in absence of cutoff. At weak coupling the 
Polyakov loop operator will remain the only massless mode and one can write down an 
effective action by integrating all other modes out. The effective action cannot be written 
in a closed form a priory, but from the requirement of gauge invariance and Gauss's law 
one can show that it will contain terms like tr(U ni )...tr(U nk )tr(U~ ni ~'"~ nk ). As evident 
from (3.17) at weak coupling mass of n-th moment of distribution function p is given by 
V n . Since at a fixed temperature V\ is the smallest, we expect that p\ = j^tr{U) would be 
the lightest mode. So at weak coupling it is possible to integrate out all other moments 
and write down an effective action in terms of p\ only. One extrapolates this to strong 
coupling and assumes so far the universal features of the system is concerned, one can 
write down an effective matrix model action quadratic in pi with coefficients depending 
on the temperature the 't Hooft coupling A. 

On the basis of the above arguments a matrix model consisting of a unitary matrix 
was introduced in [2]. The partition function and effective action are given by 

Z(X,T) = J dUe s ^f {u \ (4.1) 

S eff (U) = a(X,T)tT(U)tT(U)^ + b -^{tT(U)tT(U)^\ (4.2) 

respectively, which involves the order parameter characterising the deconfined phase of 
the gauge theory given by the expectation value of the Polyakov loop (l/N) (trU). The 
equations of motion that follows from the action can be written as, 

1 

ap + 2bp i = p < p < -, 

[ 1 < P < 1, (4.3) 



4(1 -P) 2 

where p 2 = (1 / 'N^trUtrU^ . For later use let us note that the equations of motion can 
be integrated in the form of a potential. In particular, for p > 1/2 the potential can be 
written as 

V(p) = 2ap 2 + 2bp 4 + log(l -p)+f, (4.4) 

where we have added a constant, / = log(2) — 1/2 to make the potential continuous at 
p = l/2. 

Considering this model as an effective theory in the strong coupling regime, the 
phases were studied in [2] . It was shown that this model undergoes two phase transitions 
as a function of temperature. There is a first and a third order transition when b > and 
a < 1. The thermal history matches with that of the bulk gravity as long as both a and 
b are increasing functions of temperature. The first order transition can be mapped to 
the HP transition between thermal AdS and the big black hole. It was further proposed 
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that the third order transition corresponds to the black-hole string transition for the small 
black hole. This transition is not visible in the gravity approximation of string theory. 



Since on the gravity side, turning on the cutoff does not lead to any drastic change 
in the thermodynamic behaviour it is a reasonable assumption that the above arguments 
leading to the quadratic matrix model effective action will remain valid for non-zero cutoff 
and it is possible to describe the thermodynamic behaviour with the effective action given 
in (4.2), where the coefficients a and b will depend on the value of the cutoff as well. 
In what follows, we analyse this matrix model and comparing with the results obtained 
from gravity side we study the behavior of the coefficients a and b as functions of infrared 
cutoff. For convenience we write the cutoff parameter as A IR , which is related to the 
cutoff r on the gravity side through the relation A IR l 2 ~ r . 

We now use a method similar to [7] in order to study the temperature dependence 
of a, b at various 1R cutoffs. Subsequently, we discuss, with a mild assumption, how this 
method also allow us to find A dependence of a and b at fixed temperature. The procedure 
is as follows. On the matrix model side, at a generic temperature, there are two saddle 
points arising as the solution of (4.3). Of which the unstable one (maximum) corresponds 
to the small black hole and the stable one (minimum) corresponds to the large black hole 
on the gravity side. We identify the matrix model potential at maximum and minimum 
with the gravity action evaluated with the smaller and larger values of horizon radius 
respectively. That will lead to the following two equations 

2ap\ + 2bp\ + log(l - Pl )+f= -h(r ), (4.5) 
2ap\ + 2bp\ + log(l -p 2 ) + f= -J 2 (r ), (4.6) 

where the I\${ro) are the actions for the large and small black-holes respectively and 
pi t 2 are the corresponding solutions in the matrix model. Since we need the black hole 
solutions on the gravity side and the saddle points on the matrix model these equations are 
valid at T > To. At the black hole nucleation point, where T = To, above two equations 
converges into one. Therefore we can begin from an input value of a(To) while b(To) will 
be determined by that equation. In addition to (4.5,4.6) p\ and p 2 satisfy the saddle point 
equation (4.3) as well which leads to two more equations one for each of them. 

m + 2hA = 5<rbo ' (47) 
«" + 2b »i = liT^y (4 - 8) 

This set of four equations given by (4.5,4.6,4.7,4.8) can be used numerically to eliminate 
pi, p 2 , and to solve for a and b. 2 and thus provide the temperature dependence of a and 
b. 



2 We have set I, UJ3, K5 to 1 in the numerical computations. 
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The variations of a and b as function of temperature for various values of r are 
shown in figures (5). In absence of cutoff both a(T) and b(T) are monotonically increasing 
function of temperature as shown in [2]. From figure (5) we observe that b(T) remains a 
monotonically increasing function of T for all non-zero values of the cutoff r . However, 
from figure we find that a(T) increases slowly for r < 0.1. It becomes a decreasing 
function of temperature for higher values of the cutoff. This is somewhat puzzling from 
the following point of view. At zero coupling, a(T) has an interpretation of a single particle 
partition function. This increases with temperature. Thus one expects, even at a very 
weak coupling, a(T) would remain a monotonically increasing function of temperature. 
In the absence of cutoff, we have also seen in [7] that the gravity data suggests a(T) 
continues to increase with T even in the strong coupling. In view of this one may have 
expected similar behaviour of a(T) even in the presence of a cutoff. However, what we 
find here does not match with this expectation. 
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Figure 5: Plot of a(T) and b(T) for various values of cutoff r 



In the figure (6), we plot the thermal paths in the (a — b) space with b along the 
vertical-axis. Each dot in the figure has the coordinate a(T) and b(T) for a given value 
of cutoff r . the line "N" denotes the "Nucleation" point (a (To), b(T )) of the blackholes 
in the (a — b) paremeter space. Similarly, line "D" is the "Deconfinement" transition line 
corresponds to (a(T#), 6(T#)). For a given value of cutoff, the dot moves upward with 
increasing temperature begining from the nucleation line. All dots on the nucleation line 
are at the same temperature independent of cutoff since the nucleation temperature of 
black holes in the bulk is unaffected by the cutoff while it is different on the deconfinement 
line. 

From the plot, following the dots along the path, one can see that deconfinement 
temperature increases with increasing r . From the matrix model side, it can be seen 
that as cutoff increases the deconfinement temperature will increase as long as b remains 
an increasing function. We can compare this with the zero coupling situation. There the 
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Figure 6: Thermal paths in the (a — b) space corresponding to the various values of the 
cutoff in gravity, (r ) which increases from blue towards red. The direction of flow as the 
temperature increases is from bottom towards top. The vertical axis is b. 

transition temperature is determined by a(T) = 1. As we raise the value of the cutoff, 
a(T) increases more and more slowly with temperature and as a result the transition 
temperature increases. The free case can be thought of as a flow of a(T) along the 
horizontal axis where b = 0. 

The theory in this section has three scales, the temperature (3, the radius of the 
sphere Rs3 and the cutoff R c (which is related inversely to r from the gauge/gravity 
duality), a and b are thus functions of two dimensionless quantities: (3/R c and (3/Rss. 
However it can be seen from Figure (5) that the terms containing R^ 3 /((3R C ) dominate 
over the other dimensionless quantities. For fixed -R53, both a and b scale with respect to 
(3 in the same way as the cutoff r . This dependence on two dimensionless combinations 
is in contrast to the theory without the cutoff with boundary S 3 x S 1 or the one with 
flat boundary and a cutoff. In the latter cases the free energy in gravity and hence the 
parameters a and b show much simpler dependence on the scales. 

The critical exponent around the nucleation temperature, T = Tq for this matrix 
model has been calculated in [2]. This is obtained by expanding the effective potential 
(4.4) in powers of (T — T )/T . As there is no explicit dependence of the parameters on 
the cutoff here, the critical exponents thus remain the same. This is consistent with what 
we found in gravity, that the various critical exponents are insensitive to the cutoff. 



17 



rO=0.0 
rO=0.1 
r0=0.2 
r0=0.3 
r0=0.4 



rO=0.0 + 

rO=0.1 x 

rO=0.2 * 

r0=0.3 □ 

rO=0.4 ■ 



0.2 0.4 



0.2 0.4 0.6 0.8 1 1.2 

Temperature (T) 



Figure 7: Plots of a(T) and b(T) with various values of cutoff r for the theory on R 3 x S 1 . 



We now use these above inputs to construct a model on R 3 with an infrared cutoff. 
This is the subject of the next section. 



5 Effective model at strong coupling on R 3 with IR 



Having discussed the theory on S 3 with an additional IR cutoff, let us now study its 
behavior as we take the radius of S 3 to infinity keeping the IR cutoff fixed. The effective 
theory with two parameters (a, b) serves as phenomenological model for the theory on 
S* 3 x S 1 . The radius of the S 3 acts as an infrared cutoff and there is a phase transition 
at a temperature proportional to the inverse of the radius. As the radius of the sphere is 
taken to infinity the theory is always in the deconfined phase. We have seen in section 
(2) that for a flat boundary with an IR cutoff introduced by hand also shows the same 
qualitative features as the theory on the sphere. Thus at the level of a phenomenological 
model, it is natural to expect that (4.2) should serve as good effective model for the 
theory on flat space with an IR cutoff. From the point of view of the boundary theory 
this relies on the assumption that the effective model with IR cutoff also has an expansion 
in terms of the Polyakov loop. The dimensionless quantity on which the phases of the 
theory depend is f3/R c , where R c is IR cutoff on the boundary. For the theory on the 
sphere R c = R s s. 

The main qualitative difference as seen from section (2) is that for the flat bound- 
ary, the gravity theory with cutoff shows only one black hole solution that exists for all 
temperatures, whereas the theory with the boundary S 3 has two. In the latter case, the 
euclidean version of this unstable black hole serves as the bounce for the AdS to black 
hole transition. In R 3 with cutoff, one may expect a similar phenomenon when we can 
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Figure 8: The thermal path in the (a — b) space for all values of the cutoff in gravity, 
(r ). The direction of flow as the temperature increases is from bottom towards top. The 
vertical axis is b. 

find an exact solution of the Einstein's equation with an inbuilt cutoff (instead of putting 
it by hand). The Lorenzian version of the missing bounce would be an unstable black 
hole. In the appendix, we provide such a construction for the soft wall model. 

From the matrix model side we have seen that it is impossible to write down a 
smooth potential that interpolates between two minima (corresponding to thermal AdS 
and the big black hole) without having an unstable maximum in between. Here in this 
analysis we will include this unstable maximum so that it lies in the range < p < 1/2. 
The interpretation of such a saddle point is that of a stringy solution not visible in the 
supergravity approximation. A transition to such a string phase from the black hole phase 
was interpreted in [2,9] as the gravity dual of the third order phase transition [20,21] on 
the boundary. In a complete gravity solution with cutoff we expect to see an unstable 
black hole as mentioned in the earlier paragraph the corresponding saddle point for which 
in the matrix model would lie in the region p > 1/2. 

To implement this in the following numerical analysis similar to that in section (4), 
we keep the unstable saddle point at p = 0.49. The following plots of a and b (Figure(7)) 
reflects the fact that they are functions of the dimensionless quantity (3/R c . 

Figure (8) shows the thermal path in the (a — b) space. Note that there is only one 
path for all values of as against Figure (6), where we have different paths for different 
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values of r . This is primarily due to the fact that here we have only one other scale, 
the cutoff ro or R c besides the temperature. So a change in r$ will keep us on the same 
path as it can be compensated by a change in temperature. The fact that the transition 
temperature increases with r is thus obvious as the phases depend on the dimensionless 
ratio (f3r /l 2 ). This is not so for the theory discussed in the earlier part of this section 
where apart from j3 we had both R c and R s s. In the numerical calculations here the 
critical value of (a, b) is (0.681,0.664). It lies on the deconfmement line (dashed line in 
Figure (8)). For r = this value is never reached and the thermal path of a — b always 
remains above the deconfmement line. For rest of the finite values of the path crosses 
the deconfmement line and at low enough temperatures the values of a and b correspond 
to those in the confined phase. However, since the unstable maximum is always below 0.5, 
the curve for all values of r never crosses the third-order (black hole-string) transition 
line (solid line in Figure (8)) from above. The thermal history for nonzero temperature 
thus starts beyond this line. 



6 Discussion 



The analysis in this paper rests on the assumption that near the critical temperature, the 
effective theory governing the phase transitions for M = 4 SYM on S* 3 and the same on 
R 3 with an IR cutoff are identical. Some evidence for this assumption is obtained from 
the phase transitions in gravity. On the gravitational side there is no phase transition 
for a theory with boundary, R 3 . The theory in this case is always in one of the phases; 
specifically the black hole free energy is less than that of thermal AdS for all non zero 
temperatures. However there are two ways to get a phase transition at finite temperature; 
one by changing the boundary to S 3 and the other by putting a radial cutoff in gravity 
r , that has been discussed in section (2). These result in lifting the free energy of the 
black hole above zero below a certain critical temperature. To this effect, the physical 
roles played by the cutoff r and the finite radius of the three-sphere are identical. We 
have thus proposed that the effective matrix model (4.2) should also model the phase 
transitions for the theory with the cutoff. Comparing with the thermodynamic behaviour 
of its gravity dual, we analyzed this matrix model. This model is parametrized by two 
parameters a and b. We have numerically calculated the variations of these parameters 
with respect to the temperature and the cutoff. The relative behaviors of a and b with 
respect to the scales reflect those of the free energy of gravity. While for the theory on 
S* 3 with a cutoff, due to the presence of three scales (fl,Rss and R c ), the nature of the 
relative dependences is not that clear, it is quite simple when we have only two scales {(5 
and R c ). 

Apart from the similarities of the two configurations, there are of course certain 
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differences that we find in the analysing the response of the matrix model parameters 
to changes of temperature when considering S 3 or R 3 with cutoff. For the former, the 
parameters are increasing functions of temperature (see [7]). This is not so in the latter 
case (Figure (7)). Also, from the computation of the potential for a quark-antiquark pair, 
we find that the potentials differ in the above two cases. These differences could be due 
to the details of the form of the free energy from gravity. A weak coupling computation 
with an IR cutoff on the gauge theory side may shed some light on the similarities and 
differences of the two configurations. 

Acknowledgements: We would like to thank Bala Sathiapalan for discussions and for 
his comments on the paper. We also thank Souvik Banerjee for collaboration on results 
presented in the appendix. The work of TKD is supported by the South African Research 
Chairs Initiative of the Department of Science and Technology and National Research 
Foundation. Any opinion, findings and conclusions or recommendations expressed in this 
material are those of the authors and therefore the NRF and DST do not accept any 
liability with regard thereto. 

A Dynamical cutoff? 

In (5), we saw that it was impossible to write down a smooth potential which would 
interpolate between two minima (representing thermal AdS and black hole in the presence 
of cutoff) without having an unstable maximum. This maximum, in the gravity side, 
would represent a small unstable black hole which is not seen in the hard wall model 
of [13]. Furthermore, on generic ground, one may expect that since in this model, there is 
a transition from AdS to black hole, there should be a bounce solution. The Lorentzian 
version of this bounce would be an unstable black hole. Can the absence of such a solution 
in thiese models be an artifact of introducing cutoff by hand? Within the framework of 
soft wall model, we try to address this issue in the appendix. 3 

Here, we consider a toy model in five dimensions which has a dilaton along with 
a non-trivial dilaton potential. This potential goes to a negative cosmological constant 
for small value of the scalar. We construct a family of solutions for the metric with a 
soft wall like behaviour of the dilaton. At large distance, the metric asymptotes to a flat 
AdS black hole. On the other hand, at short distance, due to dilaton backreaction, the 
metric gets deformed from it asymptotic structure. We indeed see that this deformation 
induces a small unstable black hole in the geometry. This perhaps suggests that had we 

3 This computations in this appendix was done along with Souvik Banerjee, Institute of Physics, 
Bhubaneswar, India. 
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dynamically intoduced a cutoff in the soft wall model of [13], a small black hole would 
prehaps have showed up. However, it should be noted that, we only consider a toy model 
and it requires a fine-tuned dilaton potential. 

The model we consider is represented by the following action 

S = J d b x^(^-R- i(90) 2 -y(0)). (A.l) 
We look for classical solution of this action where the metric has the following form 

dS 2 



= e 2A ( K - hde + dx i) + e — dr ^ ( A - 2 



h 

i=i 

where A, B, h are functions of r only. The equations of motion can be written down by 
varrying the action with respect to the metric and the dilaton. They are 

h" + (4A' - B')ti = 0, 

A" _ A 'B' + -(<j)') 2 = 0, 
3 

6A'ti + h(2AA' 2 - 2(f)' 2 ) + We 2B = 

✓ + > + <"'-flW-£|£ = 0. (A.3) 

Here the first three equations follow from the variation of the metric and the last one is 
the scalar equation. Differentiating the third equation in (A.3) and using the first two, 
we get the last equation. In this sense, there are only three independent equations. Our 
will look for solutions of the above equations where the scalar depends on r as 



a 2 



r 



2 



(A.4) 



and then determine A,B,h and V by solving (A.3). It is easy to see that the following 
metric satisies the equations of motion 

— 2 ° 4 2 3 

dS 2 = -r 2 (h + h ie -£)dt 2 + 6 + r 2 Y^dx 2 . (A.5) 

In the expression of the metric ho and hi are constants. The role of these constants will 
be discussed shortly. The potential V for which (A.5) is a solution is given by 



V = -2 



h 1( j) 2 e^ + (3 - <j) 2 )(h e^ + h^) . (A.6) 



We first note that for large r, as goes to zero, the potential goes to a negative constant 
given by 

^Uoo = -6(/i + /ii). (A.7) 
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Setting 

ho + hi = l/l 2 , (A.8) 

we get, asymptotically, the standard negative cosmological constant term (— 12/7 2 ) in the 
action (A.l). Furthermore, coefficient of <j) 2 in the potential, for — > 0, saturates the BF 
bound m = — 4// 2 . To leading order, the metric reduces to 



dS 2 



/r^MS dr 2 

XT 1 j i= o 



(A.9) 



where we have used relation (A.8). This is an AdS black hole solution with flat horizon 
if we identify the mass of the black hole as 



M 



hia 4 



(A.10) 



Calculating R ilvpa W vpa ', it is easy to check that (A. 5) has a curvature singularity at r = 0. 



The horizon of the metric is given by the zero of h. We denote the horizon radius 
as r+. It satisfies 

1 



r + = a 



(A.11) 



The temperature can be computed by analysing the behaviour of the Euclidean 
version of the metric (A. 5) near the horizon. For that, we define a new coordinate p as 



ho + hie — p 2 . 



(A.12) 



In terms of p, (A. 5) takes the form 



dS 2 



9r= 



Ah 2 a 8 



2 Ah 2 a s p 2 2 



dp z + 



9r 6 + 



dr' 



+ r 2 + ^dx 2 , 



(A.13) 



i=i 



where r = —it. This metric does not have a conical singularity at p = only when r has 
a periodicity 

Inverse of /? gives the temperature 



T 



(l-e 



3vr/ ; 



In the limit of large r + , 



7t/ 2 



(A.15) 



(A.16) 
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Figure 9: In this figure, we have shown the behaviour of T as a function of r + for 1 = 1, 
a = 1. Similar to the spherical AdS-Schwarzschild black hole, it shows that below a 
critical temperature, horizon does not exist. 



This is the temperature of a flat black hole in the absence of a scalar. For generic value 
of r+, the behaviour of the temperature is shown in figure (9). We first not that the 
behaviour is reminiscent to that of an AdS-Schearzschild black hole with sperical horizon. 
The temperature has a minimum when r + satisfies 

a 4 

9(e^ - l)r% - 4a 4 = 0. (A.17) 
If we denote this value of r + as r m , we find that the minimum temperature is 



T m = (A.18) 

Temperature T m is identified as the nucleation temperature where a pair of black hole 
nucleates. At any temperature above T m , there are two black holes. The smaller one is 
clearly unstable with negative specific heat - it's temperature decreases as it increases in 
size. The bigger one, on the other hand, is always stable. 
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